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ABSTRACT 



Objective: The aims of this study were to compare different surgical approaches to rapid 
canine retraction by designing and selecting the most effective method of reducing 
resistance by a three-dimensional finite element analysis. Material and Methods: Three- 
dimensional finite element models of different approaches to rapid canine retraction by 
reducing resistance and distraction were established, including maxillary teeth, periodontal 
ligament, and alveolar. The models were designed to dissect the periodontal ligament, 
root, and alveolar separately. A 1.5 N force vector was loaded bilaterally to the center 
of the crown between first molar and canine, to retract the canine distally. The value of 
total deformation was used to assess the initial displacement of the canine and molar at 
the beginning of force loading. Stress intensity and force distribution were analyzed and 
evaluated by Ansys 13.0 through comparison of equivalent (von Mises) stress and maximum 
shear stress. Results: The maximum value of total deformation with the three kinds of 
models occurred in the distal part of the canine crown and gradually reduced from the 
crown to the apex of the canine; compared with the canines in model 3 and model 1, the 
canine in model 2 had the maximum value of displacement, up to 1.9812 mm. The lowest 
equivalent (von Mises) stress and the lowest maximum shear stress were concentrated 
mainly on the distal side of the canine root in model 2. The distribution of equivalent (von 
Mises) stress and maximum shear stress on the PDL of the canine in the three models 
was highly concentrated on the distal edge of the canine cervix. Conclusions: Removal 
of the bone in the pathway of canine retraction results in low stress intensity for canine 
movement. Periodontal distraction aided by surgical undermining of the interseptal bone 
would reduce resistance and effectively accelerate the speed of canine retraction. 

Keywords: Distraction osteogenesis. Oral surgical procedures. Orthodontics. Finite element 
analysis. 



INTRODUCTION 

Currently, conventional fixed orthodontic 
treatment requires about 1-2 years 4 . More time 
is required for extraction cases, such as for adult 
patients, which is a great concern and poses high 
risk of caries 5 , external root resorption 18 . Thus, 
accelerating orthodontic tooth movement and the 
resulting shortening of the treatment duration 
would be beneficial. Many researchers have 
utilized different biochemical methods involving 
medications to improve the speed and quality of 



orthodontic treatment, but the systemic influence 
on the body's metabolism makes this difficult to 
apply in Orthodontics. Recently, investigators have 
begun studying local techniques for stimulating 
better orthodontic tooth movement. Surgically 
aided rapid tooth movement has become one of the 
novel techniques for accelerating canine retraction. 

For premolar extraction cases 25 , the first phase 
of treatment is distal movement of the canines. 
Biologic tooth movement can be achieved 23 with 
conventional orthodontic treatment techniques, 
but the canine retraction phase usually lasts for 
6 to 8 months. Extraoral or intraoral anchorage 
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mechanics are required to maintain the space 
obtained during canine distalization, particularly 
when maximum anchorage is required. To date, 
several novel surgical modalities and techniques 
for rapid canine retraction have been reported to 
accelerate canine distalization. These attempts fall 
mainly into two categories: the first, dentoalveolar 
distraction osteogenesis, evolved from distraction 
osteogenesis 28 ' 16 . Osteotomies surrounding the 
canines are created to achieve rapid movement 
of the canines in the dentoalveolar segment, 
in compliance with the principles of distraction 
osteogenesis. The second technique is periodontal 
distraction aided by surgical undermining of the 
interseptal bone 14 20 . Nevertheless, the question 
remains regarding which of the above techniques 
is more effective for rapid canine retraction. 

Until now, these research modalities and 
techniques for accelerating canine retraction have 
been applied in animal experiments and clinical 
case reports. However, the biomechanics of the 
application process has not been studied. 

Biomechanics is crucial in all aspects of 
Dentistry, especially in Orthodontics. After the 
application of an inappropriate external force to 
the oral environment, the development of a stress 
field in the supporting tissues typically leads to an 
unsuccessful outcome. Finite element analysis has 
become a powerful tool for dental biomechanical 
research due to its increased availability, capacity, 
and ease of use of computer software in biologic 
modeling. It can be considerably effective and, most 
importantly, non-invasive. 

What is the status of stress in the maxillary 
teeth, the periodontal ligament (PDL), and alveolar 
bone during rapid canine retraction? How does 
the biomechanics differ between dentoalveolar 
distraction osteogenesis and periodontal distraction 
aided by surgical undermining of the interseptal 



bone? Will it produce side-effects in addition to 
tooth movement? To address these questions via 
the finite element (FE) method, in this study we 
constructed a three-dimensional (3D) finite element 
model of the maxillary teeth, the periodontal 
ligament (PDL), and alveolar bone after extracting 
the first premolars. 

MATERIAL AND METHODS 

3D FE model generation 

For better simulation of the movement of 
canine retraction during orthodontic treatment 
with different auxiliary surgical methods, a patient 
with bimaxillary protrusion was selected from the 
orthodontic department of First Affiliated Hospital 
of Xinjiang Medical University. The patient was 
informed of the experimental protocols and 
provided informed consent. The study protocol 
and patient consents were approved by the Ethics 
Review Committee. Four maxillary first premolars 
were extracted as part of the treatment. 

After the maxillary teeth were aligned and 
leveled by means of an invisible aligner, images of 
the patient's teeth and alveolar bone were obtained 
in a DICOM (digital imaging and communication in 
medicine) data format via 64-slice (a slice thickness 
of 0.5 mm) spiral computed tomography (CT). 
The geometric shapes of the maxillary teeth and 
alveolar bone (cancellous and cortical bone) were 
reconstructed by MIMICS (version 15.0; Materialise, 
Leuven, Belgium), then exported to Geomagic 
Studio (2013, Raindrop Inc., Rock Hill, South 
Carolina, U.S.A.), and modified by Unigraphics NX 
(version 8.5, Siemens, Munich, Germany). The 
initial 3D finite element model of the maxilla and 
the maxillary teeth was established with 416,931 
nodes and 227,784 solid elements (Figure 1), and 
analyzed by ANSYS Workbench (version 13.0; 



ANSYS 
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Figure 1- The initial 3D finite element model of the maxilla and the maxillary teeth 
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ANSYS Inc., Canonsburg, Pennsylvania, U.S.A.). 

Although periodontal ligament thicknesses 
differ according to age, position, and individual 
variations, the thickness of the periodontal ligament 
was considered to be consistently 0.25 mm 3 - 21 . 
The 3D finite element models of the alveolar bone 
were fabricated to fit the teeth and the periodontal 
ligament, and the thickness of cortical bone was 
considered to be 2 mm 1317 . 

In this study, we defined the non-surgically- 
aided canine retraction model as model 1. The 
numbers of nodes and elements of the initial model 
are shown in Figure 2. According to the two types 
of auxiliary surgical methods for accelerating canine 



retraction, we simulated the same surgical process 
in the initial 3D finite element model, and two 3D 
FEMs were generated based on the two different 
surgeries for comparison with the initial model 
(model 1). 

We defined the FEM of periodontal distraction 
aided by surgical undermining of the interseptal 
bone as model 2 (PD). After the first premolar 
was removed, the interseptal bone mesial to the 
extraction socket was undermined. The bone in the 
premolar socket was eliminated vertically along 
the buccal and lingual sides, extending obliquely 
toward the base of the interseptal bone to weaken 
resistance. The interseptal bone was not cut 



Models 


Number of nodes 


Number of elements 


Cancellous bone 


69.897 


41.814 


Cortical bone 


118.762 


68.648 


Canine 


7.293 


4.181 


Second premolar 


3.933 


2.241 


First molar 


5.961 


3.497 


PDLof canine 


17.383 


8.584 


PDL of second premolar 


11.482 


5.645 


PDLof first molar 


18.873 


9.306 



Figure 2- Number of nodes and elements generated for initial model 
PDL=periodontal ligament 




Oblique undermining groove 



Vertical undermining groove 




Figure 3- Surgical techniques for undermining interseptal bone distal to the canine. No cuts are performed on buccal and 
lingual plates. Note depths and positions of undermining grooves 




Removal of buccal cortical 
bone of extraction socket 



Outlined canine root with 
the dentoalveolarsegment 




Figure 4- Surgical techniques for dentoalveolar distraction osteogenesis 
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through mesio-distally toward the canine (Figure 
3). The depth of the undermining grooves was 
dependent on the thickness of the interseptal bone. 

The final model (model 3) was the FEM model of 
dentoalveolar distraction osteogenesis (DDO). The 
first premolar was extracted, and the buccal cortical 
bone of the extraction socket was removed. An 
osteotomy line was determined between the buccal 
root apex of the first premolar and canine apex 
(Figure 4). The root of the upper canine tooth was 
outlined mesially and distally, with the dentoalveolar 
segment to be used as a transport disc at the apical 
region. The cortical bone, especially at the level of 



the apical region, was also eliminated for maximal 
bodily movement during distraction. 

Material properties 

The mechanical properties of the teeth, PDL, 
and cortical and cancellous bone were defined 
as isotropic, homogenous, and linearly elastic. 
Young's modulus and Poisson's ratio for all materials 
were calculated according to previously described 
methods (Figure 5) 19,24 . The boundary condition 
was defined as the exterior border of the maxilla. 



Material 


Young's Modulus (MPa) 


Poisson's Ratio 


Cancellous bone 


1.370 


0.3 


Cortical bone 


13.700 


0.3 


Tooth 


20.000 


0.3 


PDL 


0.68 


0.3 



Figure 5- Mechanical properties of the materials used in the finite element models 
PDL=periodontal ligament 




Figure 6- The coordinate system. A mesiodistal simulated retraction force of 1 .5 N was loaded bilaterally to the center of 
the crown between the first molar and the canine 
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Experimental conditions 

The coordinate system and directions of forces 
and loads are shown in Figure 6. A simulated 
retraction force of 1.5 N (approximately 150 g) 
was loaded bilaterally to the center of the crown 
between the first molar and the canine, and stress 
distribution and its magnitude were analyzed by 
ANSYS Workbench, a 3D finite element analysis 
program. An assessment of the stress on the bone 
elements was performed by von Mises equivalent 
stress and maximum shear stress, and we used 
the total deformation value to assess the initial 
displacement of the teeth and the PDL. 



Table 1- Total Deformation values induced in three models 



RESULTS 

Canine 

Immediately after loading of the retraction 
forces of 1.5 N (approximately 150 g) bilaterally 
between the crown of the first molar and the canine, 
the initial displacement of the canine was highly 
concentrated in the distal area of the crown in all 
three models. Compared with the total deformation 
of the canine in model 1 (non-surgical) and model 
3 (DDO), the canine in model 2 had the maximum 
value of total deformation of 1.9812 mm (Table 1, 
Figure 7). 

Compared with the maximum value of equivalent 



Total Deformation values (mm) 




Canine 


First molar 


PDL of Canine 


PDL of First molar 


Non-surgical 


1.8925 


0.15441 


0.10237 


0.070115 


PD 


1.9812 


0.12674 


0.13265 


0.069375 


DDO 


1 .9426 


0.14076 


0.11349 


0.069725 



PD=periodontal distraction; PDL= periodontal ligament 
DDO=dentoalveolar distraction osteogenesis 



2.5 




I Non-surgical 
IPD 
DDO 



Canine First molar PDL of Canine PDL of First 

molar 

Figure 7- Total deformation values induced in different parts of three models. PDL= periodontal ligament 
Table 2- Equivalent (von Mises) Stress values induced in three models 



Maximum Von Mises stress values (MPa) 





Canine 


First molar 


PDL of Canine 


PDL of First 
molar 


Cortical bone 


Non-surgical 


0.56251 


0.57389 


0.049053 


0.035727 


1,7942 


PD 


0.35494 


0.43203 


0.043629 


0.03509 


0.622 


DDO 


0.45463 


0.45208 


0.043347 


0.035098 


0.6796 



PD=periodontal distraction; PDL= periodontal ligament 
DDO=dentoalveolar distraction osteogenesis 
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(von Mises) stress on the canine in models 1 and 3 
(DDO), the canine in model 2 (PD) had the lowest 
equivalent (von Mises) stress value of 0.35494 MPa 
(Table 2, Figure 8). Likewise, compared with model 
1 (non-surgical) and model 3 (DDO), the value of 
maximum shear stress of the canine in model 2 
(PD) was the lowest (Table 3). 



Molar 

As the anchor tooth, the value of total deformation 
on the first molar in model 2 was 0.12674 mm, 
which is the lowest in comparison with that of the 
molar in model 1 (non-surgical) and model 3 (DDO) 
(Table 1). 

Distribution of maximum equivalent (von Mises) 
stress on the root and in the PDL of the first molar 



Table 3- Maximum Shear Stress values induced in three models 



Maximum Shear Stress values (MPa) 





Canine 


First molar 


PDL of Canine 


PDL of First 
molar 


Cortical bone 


Non-surgical 


0.28301 


0.29064 


0.028057 


0.020853 


2,0781 


PD 


0.23029 


0.22253 


0.025144 


0.020169 


0.37052 


DDO 


0.25017 


0.25022 


0.024981 


0.020165 


0.62196 



PD=periodontal distraction; PDL= periodontal ligament 
DDO=dentoalveolar distraction osteogenesis 




i Non-surgical 
IPD 
DDO 



Canine First molar PDL of PDL of First Cortical 
Canine molar bone 

Figure 8- Equivalent (von Mises) stress values induced in different parts of three models. PDL= periodontal ligament 
2.5 



1.5 



in 

E 
3 

E 



0.5 




i Non-surgical 
IPD 
DDO 



Canine First molar PDL of PDL of First Cortical 
Canine molar bone 



Figure 9- Maximum shear stress values induced in different parts of three models. PDL= periodontal ligament 
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in all three models was concentrated in the root 
furcation of the molar, and the value of maximum 
equivalent (von Mises) stress of the first molar in 
model 2 (PD) was the lowest (Table 2). The value 
and distribution of maximum shear stress on the 
root and in the PDL of the first molar in all three 
models showed the same tendency as the maximum 
equivalent (von Mises) stress, although the values 
differed. 

Alveolar 

Compared with cancellous bone, the buccal side 
of cortical bone around the canine root was the most 
stressed area. The value of maximum equivalent 
(von Mises) stress on the cortical bone around the 
canine root in model 2 (PD) and model 3 (DDO) was 
lower than that in model 1 (non-surgical) (Table 2). 
The value of maximum shear stress on the buccal 
side of cortical bone around the canine root in model 
2 (PD) was the lowest (Table 3). 

DISCUSSION 

In this study, orthodontic force was applied to 
three FEMs which simulated two different surgical 
interventions and conventional treatment for canine 
retraction. Stress distribution and deformation on 
the root, PDL, and cortical bone were evaluated. 
The stress patterns in the PDL, root, and bone 
were displayed in separate illustrations, so that 
the complex responses of these tissues to different 
types of surgical interventions could be readily 
compared. The forces were applied to the surfaces 
of the teeth, mesio-distally, as in normal clinical 
practice. The results of this study, which showed 
stress distribution along the root, PDL, and alveolar 
bone, provided insight into clinical observations. 

Orthodontic tooth movement is a biological 
process characterized by sequential reactions of 
periodontal tissue against a biomechanical force 
system 6 . It is also a process in which the application 
of a mechanical force induces alveolar bone 
resorption on the pressure side and alveolar bone 
deposition on the tension side 11 . The orthodontic 
force system is a complicated three-dimensional 
system which is difficult to evaluate in clinical 
conditions and orthodontic force plays an important 
role in the entire biomechanical process during 
tooth movement. In our study, it was found that the 
effects of force application would be changed when 
the force-loading environment was changed. For 
accelerating canine retraction speed, periodontal 
tissue, especially hard tissue such as alveolar bone 
around the canine, is the most important source of 
resistance. 

Compared with the conventional method of 
canine retraction, we found that DDO (dentoalveolar 
distraction osteogenesis) and PD (periodontal 



distraction aided by surgical undermining of the 
interseptal bone) could reduce resistance during 
canine retraction in this finite element analysis. This 
finding was consistent with those of previous clinical 
reports 10 - 14 and animal experiments 7 ' 20 with rapid 
orthodontic tooth movement. Compared with DDO, 
it was noteworthy, from the mechanical analysis in 
this study, that PD had more advantages in reducing 
resistance and accelerating canine movement. It 
also showed that the stress along the distal side 
of the canine root was reduced more significantly 
in model 2 (PD) (Figures 8, 9), with more obvious 
canine displacement than in model 3 (DDO) (Figure 
7). The main pathway of canine movement was 
via the cancellous bone and the interseptal bone, 
since the main portion of cancellous bone in the 
pathway of canine distalization was removed by PD. 
As a result, it is more effective to accelerate canine 
retraction through periodontal distraction aided by 
surgical undermining of the interseptal bone. 

On the basis of previous studies, resistance to 
tooth movement is increased when the roots are 
torqued lingually or buccally. This principle was used 
by Rickets 22 (1979) and is called cortical anchorage. 
Generally, cortical bone offers more resistance to 
resorption. The cortical bone could also block tooth 
movement in most cases in orthodontic treatment. 
In this study, we observed the distribution of stress 
on the buccal side of cortical bone around the canine 
root, which implied that the cortical bone on the 
buccal side of the canine was also the source of 
resistance to canine movement. That explained the 
principle of surgical procedures for DDO, which was 
designed with the dentoalveolus as a bone transport 
segment for posterior movement 10 . However, based 
on our results, PD was more effective in canine 
retraction. Furthermore, in clinical applications, the 
interseptal bone distal to the canine was undermined 
to weaken its resistance immediately after dental 
extraction, and a tooth-borne appliance was used 
for rapid tooth retraction. This procedure would 
minimize surgical injury and reduce complications 
for patients. 

In addition to proper orthodontic force, the 
center of resistance (CR) of a tooth is also a critical 
factor in predicting and planning the esthetic 
movement of anterior teeth 12 . Ideally, force should 
be applied as closely as possible to the CR of a tooth, 
to achieve movement of the tooth in orthodontic 
treatment. Although we cannot apply the force on 
the CR, we could change the conditions around 
the CR. The concentration point of stress on the 
pressure side of the canine root tended to move 
closely to the crown after surgical undermining of 
the interseptal bone in the PD model. In that case, 
the canine CR became closer to the crown through 
resistance-reducing surgery, especially in the PD 
model, and this would lead to possible canine 
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body movement. Otherwise, the lowest maximum 
equivalent (von Mises) stress value on the distal 
side of the root surface in model 2 (PD) implied 
that reducing the resistance around the canine CR 
in that model would reduce the stress value. The 
lower and more uniform stress on the root indicated 
that translational tooth movement may be achieved, 
which made it more effective in accelerating canine 
retraction speed. 

However, it should be pointed out that the 
concentrated point of force loading at the distal side 
of the canine cervix for the PDL and alveolar bone 
shifted to the lingual side in the PD model, and, 
conversely, in the DDO model, this point shifted to 
the buccal side. This indicated to us that the canine 
had a tendency to rotate lingually in the PD model 
and buccally in the DDO model after resistance was 
reduced. The reason for this might be that after 
the different positions of resistance source were 
reduced by surgery, the canine would move to the 
position with the lowest resistance. Therefore, we 
should avoid this side-effect when using PD or DDO 
during orthodontic treatment. 

Anchorage control should be considered when 
accelerated tooth movement is necessary. An 
obvious strategy for anchorage control would be 
to concentrate the force needed to produce tooth 
movement where it was desired, and then to 
dissipate the retraction force in the PDL of anchor 
teeth as much as possible. In this study, the value 
of maximum equivalent (von Mises) stress in the 
PDL of the first molar in two models was surgical 
reduction of resistance to a level lower than in the 
canine. Further, the value of total deformation of the 
first molar in the DDO and PD models was far below 
that of the canine at the initial stage of force loading. 
This indicated that reducing resistance by DDO or 
PD to accelerate canine retraction was a safe way to 
protect anchorage and would not reduce anchorage 
during canine retraction. Furthermore, the pattern 
of force distribution in the PDL and root of the first 
molar showed that high stress concentration is 
observed on the root surface at the furcation level, 
which is in contrast to the canine, which displayed 
a greater apical concentration. This may explain 
the high incidence of apical root resorption 25 that 
appears to occur in the maxillary molar on the root 
surface at the furcation level, which is difficult to 
see radiographically in orthodontic treatment. 

This study and others have demonstrated that 
FEM provides a solid, workable foundation for 
modeling a system of orthodontic tooth movement 1 . 
The chief advantage of FEM is that it can be 
magnified nearly infinitely, in terms of both the actual 
volumetric construction itself and the mathematical 
variability of its material parameters. However, as 
with any theoretical model of a biological system, 
there are limitations. The mechanical behavior 



of the materials was assumed to be linear elastic 
(homogeneous and isotropic), and the value of 
each material was inferred from previous reports. 
Cortical bone thickness and cancellous bone quality 
were not incorporated into the analysis, to prevent 
bone stress from being dominated by bone quality 
and potentially confounding the outcomes related 
to other relevant factors. In addition, the stress 
analysis of soft tissues was not considered in 
this study. The soft tissues, such as gingival and 
facial muscles, are also sources of resistance for 
blocking rapid tooth movement. Regardless of these 
limitations, we integrated a finite element approach 
with variable analysis to investigate the comparative 
influences of resistance source, the pathway of 
canine movement, and different types of surgeries 
for rapid canine retraction by reducing resistance. 

CONCLUSIONS 

In this study, compared with dento-alveolar 
distraction osteogenesis, periodontal distraction 
aided by surgical undermining of the interseptal 
bone would reduce the resistance in the pathway 
of canine movement more effectively to accelerate 
canine retraction speed in rapid canine movement 
during orthodontic treatment. 

The results indicated that rapid canine retraction 
aided by the surgical reduction of resistance might 
create side effect: 1) it might lead to canine rotation 
during distalization, 2) the geometry of the root of 
maxillary first molar makes it less prone to apical 
stress concentration, therefore, resorption on the 
root surface at the furcation level of the molars is 
more easily caused. 
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